数学 启发 算法 求解 单 源 设施 区 位 及 其 变种 问题 
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摘要 : 考虑 设施 服务 区 空间 连续 和 设施 数量 限制 ， 对 经 典 的 单 源 设施 区 位 问题 
(SSCFLP) 进 行 扩展 ， 形 成 变种 问题 : 设施 服务 区 空间 连续 SSCFLP(CFLSAP)、 设 施 数量 
限制 的 SSCFLP (SSCKFLP) 和 设施 服务 区 空间 连续 及 数量 限制 的 SSCFLP(CKFLSAP)。 
针对 SSCFLP 及 其 变种 问题 ， 提 出 了 一 个 数学 启发 算法 。 该 算法 从 一 个 初始 解 开始 ， 和 迭 
代 地 对 当前 解 进行 超大 邻 域 搜 索 改 进 ， 直 到 若干 次 尝试 不 能 该 改进 当前 解 为 止 。 超 大 邻 
域 定义 步骤 如 下 : 随机 选择 一 个 客户 ， 从 当前 解 中 选择 客户 附近 O 个 最 邻近 设施 及 其 客 
户 ， 再 挑选 这 些 客户 的 最 邻近 候选 设施 。 使 用 所 选择 候选 设施 子 集 和 客户 子 集 ， 构 建 子 
问题 数学 模型 , 求解 模型 , 并 使 用 模型 解 更 新 当前 解 。 构 造 2 组 案例 数据 测试 本 文 算法 ， 
结果 表明 : 数学 启发 算法 能 够 有 效 地 求解 SSCFLP 及 其 变种 问题 ， 求 解 结果 极其 接近 问 
题 最 优 解 或 目标 值 下 界 , 相对 差异 为 0.01%(SSCFLP)、0.22%(CFLSAP)、0.00%(SSCKFLP) 
和 0.08%(CKFLSAP). 4b, SSCFLP 或 SSCKFLP 增加 设施 服务 区 空间 连续 约束 后 ， 目 
标 值 增加 幅度 不 大 ， 但 最 优 设施 区 位 可 能 发 生变 化 。 
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1 引 


Il 


设施 区 位 问题 (FLP) 是 一 大 类 寻找 最 优 设施 位 置 的 问题 。 该 问题 广泛 应 用 于 学 校 、 医 
疗 、 人 公安、 养老、 应急、 物流 等 公共 或 商业 设施 选 址 决策 。 根 据 应 用 场景 的 不 同 ，FLP 
类 型 众多 : 连续 空间 区 位 或 离散 区 位 、 有 无 设施 容量 限制 、 需 求 是 否 允 许 拆 分 指派 、 是 
否 考虑 设施 成 本 、 是 否 指定 设施 数量 等 ， 规 划 目 标 划分 为 效率 目标 、 公 平 目 标 等 。 根 据 
设施 服务 与 需求 是 否 存在 不 确定 性 ， 又 可 将 FLP 划分 为 确定 性 FLP 和 随机 性 FLP。 其 
中 ， 单 源 CFLP 问题 (SSCFLP) 是 需求 不 允许 拆 分 指派 的 、 有 容量 约束 的 、 顾 及 设施 成 本 
的 离散 区 位 问题 ， 求 解难 度 较 高 ， 也 是 近年 关注 度 最 高 的 区 位 问题 之 一 。 

令 集合 /为 潜在 设施 区 位 ,集合 ./ 为 客户 区 位 ,设施 i(i ED 有 最 大 服务 容量 s; 和 固定 建 
设 成 本 为 f， 客 户 (jE 有 服务 需求 4;， 被 设施 i 服务 的 费用 为 cy。 定义 布尔 型 决策 变量 y; 
为 是 否 在 区 位 ;建设 施 , 布尔 型 决策 变量 wj 为 是 否 指 派 设施 i 服务 客户 j，SSCFLP 数 学 模型 
如 下 : 


Minimize ie fi + Dier X jej CijXij (1) 
Subject to: Vier xij = 1,vi EJ (2) 
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È jez djXij = sip Vİ EI (3) 
yi = {0,1}, viel (4) 
xij = {0,1}, vie Lj EF (5) 
模型 中 , 目标 函数 (1) 最 小 化 设施 固定 成 本 和 获取 服务 成 本 。 约束 (2) 确 保 每 个 客户 由 唯 
的 设施 提供 服务 ; 约束 (3) 要 求 指派 给 每 个 设施 的 客户 需求 量 不 超过 其 最 大 服务 容量 。 约 
束 (4) 和 (5) 定 义 两 类 布尔 型 决策 变量 。 

求解 SSCFLP 问题 的 方法 划分 为 两 大 类 : 精确 算法 和 启发 式 算 法 (Basu et al. 2015; 
Ulukan & Demircioğlu 2015)。 精 确 算法 包括 : 分 枝 定 界 (Neebe & Rao 1983; Holmberg et 
al. 1999)、 列 生成 (Diaz & Fernandez 2002)、Cut-and-Solve (Yang et al. 2012; Gadegaard et 
al. 2018) 等 。 启 发 式 算法 包括 : 拉 格 朗 日 松弛 启发 LH)(Barcelo & Casanova 1984; 
Klincewicz & Luss 1986; Beasley 1993; Sridharan 1993; Agar & Salhi 1998; Hindi & Pienkosz 
1999; Cortinhal & Captivo 2003)、 禁 忌 搜索 ( Filho & Galvao1998; Delmaire et al. 1999; 
Cortinhal & Captivo 2003)、 超大 规模 邻 域 搜 索 (VLNS) (Ahuja et al. 2004; Tran et al. 2017)、 
4y 48 ZB (Contreras & Diaz 2008 )、 核 搜索 (Kernel Search)(Guastaroba & Speranza 2014), 
“ 廊 道 "方法 (Corridor Method)(Caserta & Voß 2020)、 松 弛 自 适 应 记忆 规划 (RAMP)(Oliveira 
et al. 2020) 等 。 因 SSCFLP 问题 具有 极 高 的 计算 复杂 度 ， 精 确 算 法 难以 高 效 求解 规模 较 
大 的 案例 ; LH 算法 速度 快 ， 但 求解 质量 不 能 保证 ; 超大 规模 邻 域 搜索 、 核 搜索 和 “万 道 ” 
方法 等 元 启发 式 算法 在 求解 质量 方面 具有 竞争 力 ， 但 计算 时 间 仍 然 偏 长 。 

在 SSCFLP 应 用 实践 中 ， 也 存在 需求 特定 需求 ， 如 拟 建 设施 数量 约束 (Aardal et al. 
2015; Wang 2017)、 要 求 每 个 设施 服务 区 空间 连续 等 。 例 如 ， 某 一 区 域 进 行 义务 教育 学 
校 布 局 规划 ， 需 要 考虑 学 校 建 设 成 本 、 学 校 数 量 要 求 ， 在 确定 学 校区 位 的 同时 划 定 其 学 
区 。 某 一 城市 进行 基层 卫生 服务 中 心 布局 规划 ， 遵 照 国家 行业 规范 ， 每 个 中 心 履 盖 一 定 
的 空间 范围 ， 满 足 15 分 钟 出 行 要 求 ， 且 考虑 中 心 的 服务 能 力 。 为 求解 这 些 问 题 ， 有 必 
要 扩展 SSCFLP 模型 ， 增 加 这 两 个 约束 条 件 。 在 分 区 问题 研究 中 ， 空 间 连 续 约 束 有 三 个 
建 模 方法 : 生成 树 模型 、 次 序 模型 和 流 模型 (Shirabe 2009; Duque et al. 2011)。 树 模型 中 
约束 条 件数 量 庞大 ， 仅 适用 于 小 规模 案例 ( 结 点 数量 小 于 50); 次 序 模型 和 流 模型 的 约束 
条 件 相对 较 少 。 针 对 p-regions 问题 ， 三 个 模型 在 3 小 时 内 最 多 能 够 求解 49 个 空间 单元 
的 案例 。 为 SSCFLP 扩展 空间 连续 约束 ， 精 确 方法 的 求解 能 力 仍 是 未 知 数 。 

基于 以 上 讨论 ， 本 文 扩展 SSCFLP 模型 ， 考 虑 设施 服务 区 的 空间 连续 性 以 及 设施 数 
量 要 求 ， 设 计 一 个 数学 启发 算法 求解 SSCFLP 及 其 扩展 问题 ， 并 测试 算法 的 求解 性 能 。 
2 问题 定义 

一 个 特定 的 地 理 区 域 ， 包 含 半 个 空间 单元 。 每 一 个 空间 单元 7 具有 属性 dj, RRE 
间 单 元 j 的 服务 需求 数量 ; n 个 单元 中 有 m(m<n) 个 单元 适合 建设 设施 , 作为 候选 设施 区 
位 ,设施 i 服务 容量 为 s;， 固 定 建 设 成 本 为 fi。 单 元 i 被 设施 j 服 务 的 费用 为 变量 cy( 空 间 
单元 i 与 j 之 间 的 单位 交通 成 本 与 需求 量 4 的 乘积 )。 令 集合 IT 表示 m 个 候选 设施 单元 ， 
集合 /为 n 个 需求 单元 (1 CJ), I SSCFLP 模型 (1)-(5) 也 适用 于 该 区 域内 设施 选 址 。 


要 求 设施 服务 区 空间 连续 ， 需 要 在 SSCFLP 模型 中 增加 空间 连续 约束 。 将 地 理 空间 
表达 为 一 个 网 络 图 ， 即 将 空间 单元 作为 结 点 ， 相 邻 空 间 单元 相连 接 作为 网 络 边 ， 形 成 一 
个 无 向 网 络 图 。 基于 网 络 图 , 可 采用 “网 络 流 ” 概 念 表 达 分 区 空间 连续 性 (Duque et al. 2011)。 
在 一 个 设施 服务 区 内 ， 每 个 单元 产生 一 个 单位 的 流量 ， 通 过 网 络 边 流动 ， 最 终 汇 入 到 区 
内 设施 单元 。 设 施 单元 不 产生 流量 ， 但 汇 入 最 高 达 n-K 单位 的 流量 (K 为 设施 数量 )。 产 
流 和 汇流 单元 所 构成 的 区 域 是 一 个 空间 连续 的 区 域 。 令 集合 Nj 表示 单元 j 的 相 邻 单元 ， 
整 型 变量 fi 表示 设施 i 服务 区 内 单元 j 到 单元 的 流量 ， 则 可 构建 空间 连续 约束 条 件 : 


fix Sn*xiviElje]J,keN; (6) 
fijx SN* XR, ViE Lj EJ,K EN; (7) 
Xren; fijk — Leen; fikj 2xipViEljEes\i (8) 
fir 2 0,ViELjEJ,k EN; (9) 


ANE (6) FT) HARE F EARE K AEEA BEE AF hE i 不 
服务 单元 j(xi; = 0)， 则 单元 j 在 i 服务 区 中 不 会 产生 流 fijx=0， 也 不 会 有 汇流 fixj=0; 当 
设施 i 服务 单元 j 在 (xi; = 世 时 ， 可 能 会 有 进出 单元 7 的 流 。 若 单元 7 不 是 汇集 单元 (设施 
单元 ), 约束 (8) 保 证 该 单元 至 少 产生 一 个 单位 的 流 。 约束 (6)-(8) 将 使 每 一 个 设施 服务 区 中 
的 产 流 汇集 到 设施 单元 ,一 个 设施 服务 区 内 , 每 个 空间 单元 必需 产生 流量 流向 邻 域 单元 ， 
最 终 汇 集 到 设施 单元 ， 这 使 得 一 个 设施 服务 区 内 任意 两 个 单元 都 是 连通 的 ， 从 而 保证 了 
服务 区 空间 连续 。 
需要 指定 设施 数量 时 ， 增 加 条 件 (10)， 指 定 设施 数量 为 Ko 

dier Vi = K (10) 
基于 以 上 约束 ,有 三 个 SSCFLP 扩展 模型 : 四 设施 服务 区 空间 连续 约束 的 SSCFLP， 
记 为 CFLSAP， 模 型 包括 式 (1)-(9); 包 设 施 数量 约束 的 SSCFLP， 记 为 SSCKFLP， 模 型 
包括 式 (1)-(5) 和 (10); 四 设施 服务 区 空间 连续 约束 和 设施 数量 约束 的 SSCFLP， 记 为 
CKFLSAP， 模 型 包括 式 (1)-(10)。 


2 数学 启发 算法 
局 部 搜索 算法 使 用 邻 域 搜索 用 于 渐进 改进 当前 解 ， 广 泛 应 用 于 离散 优化 问题 求解 。 
区 位 问题 中 ， 设 施 的 变动 必然 会 引起 需求 指派 变化 ， 使 得 邻 域 搜索 较为 复杂 。 因 此 ， 区 
位 问题 算法 设计 的 一 个 核心 问题 是 : 如 何 有 效 地 移动 设施 并 更 新 指派 ， 寻 找到 可 行 的 、 
更 节约 的 邻 域 解 。 因 设施 容量 约束 ， 涉 及 较 多 设施 的 超大 规模 邻 域 搜索 才 有 较 高 可 能 性 
改进 当前 解 。 然 而 ，SSCFLP 的 超大 规模 邻 域 结 构 定 义 复杂 ， 邻 域 搜 索 的 计算 复杂 度 也 
很 高 ， 存 在 指数 级 数量 的 客户 移动 移动 和 设施 移动 (Ahuja et al. 2004; Tran et al. 2017)。 
本 文 拟 采用 数学 模型 求解 超大 邻 域 。 算 法 原理 如 下 : OEH LH 方法 (Holmberg et al. 
1999)、 线性 松弛 或 启发 式 算法 获得 问题 初始 解 ; 包 从 当前 解 中 随机 选择 一 个 超大 邻 域 ， 
获得 邻 域内 需求 点 、 当 前 设施 和 候选 设施 集合 ; @@ 构 造 邻 域 内 子 问题 模型 并 求解 ， 更 新 


当前 解 ; ORRPVT PROMS, AIET AA EE AH. PROE PK 
邻 域 ， 步 又 @) 用 于 发 现 邻 域内 的 最 优 解 。 基 于 这 一 思路 ， 数 学 启发 算法 示意 如 下 : 
参数 :连续 未 更 新 最 好 解 循环 数 (mloops) 


1. s= GenerateInitialSolution(); 


2. notImpr=0; 

3. While notImpr < mloops: 

4. Select a large neighborhood /* and J* from solution s; 
5 s*= SolveSubProblem(/*, J*); 

6. s ’=UpdateSolution(s,s*); 

Ts s’’=VNDSearch(s’) 

8 If fis’’)<f(s): s=s’’, notImpr=0; 

9 else: notiImpr+=1; 

10. Output(s). 

以 上 算法 中 ， 步 又 (1) 获 得 初始 解 。 因 很 多 案例 难以 获得 可 行 的 初始 解 ， 本 文 按 以 下 
顺序 考虑 解 的 可 行 性 : 指派 约束 、 设 施 数量 约束 、 空 间 连 续 约 束 和 设施 容量 约束 。 考 虑 
到 设施 容量 约束 难以 满足 ， 本 文 算 法 允许 设施 容量 超 限 ， 但 在 目标 函数 中 对 超 限量 进行 
惩罚 。 增 加 整 型 决策 变量 H; = {0,1,2...}, vi ET， 表 示 超 出 设施 ; 服务 能 力 的 量 ， 容 量 约 
束 和 目标 函数 可 改 为 : 


2jeJaixij < siyi + Hy Vi EI (11) 
Minimize: ier fi + Lier Yjer CijXij + a Lier Hi (12) 
通过 以 上 模型 修改 ， 允 许 设施 容量 超 限 ， 容 易 产生 初始 解 。 设 置 设施 容量 超 限 惩罚 系数 
AWK, 后 续 算 法 逐步 减 小 Hi 值 , 直到 变量 值 为 零 , 从 而 使 容量 约束 条 件 (3) 得 到 满足 。 
这 一 设施 容量 软 约束 也 使 MIP 优化 器 快速 获得 可 行 解 ， 有 利于 模型 精确 求解 。 
初始 解 生成 方法 可 选择 : 拉 格 朗 日 松弛 启发 、 线 性 松弛 启发 、ADD 或 DROP 启发 
方法 。 拉 格 朗 日 松弛 启发 算法 (Holmberg etal. 1999) 可 直接 求解 SSCFLP 案例 ， 通 过 调整 
设施 数量 可 用 于 求解 SSCKFLP。 线 性 松弛 启发 是 将 修改 布尔 变量 xij; 为 [0,1] 连 续 变 量 ， 
使 用 优化 器 对 模型 求解 ， 对 求解 结果 进行 修复 。ADD 启发 或 DROP 启发 是 经 典 的 构造 
算法 , 以 贪心 方式 获得 可 行 解 考虑 服务 区 连续 时 , 上 述 算法 结果 需要 进行 连续 性 修复 。 
算法 步骤 (4) 从 当前 解 中 随机 选择 一 个 超大 邻 域 ， 包 括 设施 子 集 天 和 需求 子 集 J*。 
选择 邻 域 是 算法 设计 的 一 个 关键 点 , 步骤 如 下 : 也 随机 选择 一 个 需求 点 , 以 该 点 为 中 心 ， 
从 当前 解 设 施 集合 中 自 近 而 远 地 选 择 O Ai ORRE O 个 设施 的 需求 点 构成 需求 
ATEJ” OTR .7* 中 每 个 需求 点 最 近 的 设施 以 及 已 选择 的 2 个 设施 构成 子 集 I. A 
上 步骤 所 选择 邻 域 的 大 小 受 2 值 影响 ，O 值 取 值 范围 为 [Oo Omaxl> FE. Onin = 
min(L/2,7), Qmax =min(L,10), L 为 当前 解 中 的 设施 数量 。 当 Z<10 IN, O 处 于 区 间 [Z/2， 
L]; “4 11SLS13 时 ，O 处 于 区 间 [Z2,10]; “4 L>14 时 ，O 处 于 区 间 [7,10]。 
邻 域 规模 影响 算法 性 能 : 大 邻 域 大 小， 当前 设施 区 位 往往 难以 更 新 ; 若 邻 域 过 大 ， 
模型 计算 效率 低 。 当 集合 的 规模 偏 大 时 ， 可 能 会 使 模型 计算 效率 低下 。 为 此 ， 需 要 精 
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简 一 部 分 候选 设施 ， 再 进行 子 问 题 建 模 和 求解 。 精 简 方法 : 从 1* 选 择 O 个 当前 解 中 已 使 
用 设施 ， 再 从 其 余 设 施 中 随机 选择 O 个 设施 ， 更 新 集合 IY， 使 其 包括 20 个 设施 。 

步 又 (5) 使 用 设施 子 集 1* 和 需求 子 集 .1* 构 造 子 问题 模型 ， 并 使 用 MIP 优化 器 求解 。 
针对 SSCFLP、CFLSAP、SSCKFLP 和 CKFLSAP 四 个 问题 , 前 两 者 构造 SSCFLP 模型 ， 
后 两 者 构造 SSCKFLP 模型 ， 保 证 条 件 (10) 得 到 满足 。 

步骤 (6) 使 用 局 部 解 s* 更 新 当前 解 ， 获 得 新 解 i"。 对 于 CFLSAP 和 CKFLSAP 问题 ， 
新 解 s "可 能 不 满足 空间 连续 要 求 ， 此 时 需要 进行 对 新 解 进 行 修复 。 

步 又 (6) 修 复 过 程 通常 会 使 目标 值 较 差 , 步骤 (7) 使 用 可 变 邻 域 下 降 (VND) 算 法 改进 当 
前 解 。VND 算法 中 采用 1 单元 移动 和 2 单元 移动 算 子 (Kong et al. 2017)， 获 得 局 部 最 优 
解 。 

步 又 (8) 尝 试 更 新 当前 解 。 若 新 解 s”" 优 于 当前 解 s， 更 新 当前 解 ; 和 否则， 维持 当前 解 
不 变 。 

迭代 执行 步骤 (4)-(9)， 当 连续 若干 次 循环 未 能 更 新 当前 解 ， 算 法 终止。 


3 算法 测试 


3.1 案例 设计 


本 文 基于 两 个 典型 区 域 构造 测试 案例 。 区 域 ZY 为 城区 , 约 13.4 平 方 干 米 , 包含 326 
个 空间 单元 ， 有 15 所 小 学 、3783 名 学 生 ; KR ZY 为 县 域 ， 约 1000 FARK, AF 
1276 个 空间 单元 ， 有 20 个 乡镇 街道 办 ， 人 口 819812 A. Kit ZY 将 学 生 人 数 作为 需求 
量 ， 在 现 有 学 校 基 础 上 手工 增加 部 分 点 位 作为 潜在 设施 区 位 ， 设 施 容 量 参照 现 有 学 校 均 
值 。 区 域 GY 将 人 口 数量 作为 需求 量 ， 在 现 有 乡镇 街道 办 所 在 地 基础 上 手工 增加 部 分 区 
位 作为 潜在 设施 区 位 ， 设 施 容 量 参 考 现 有 乡镇 街道 办 人 口 数 量 均 值 。 图 1 为 案例 区 域 服 
务 需 求 和 候选 服务 设施 的 空间 分 布 示意 图 ， 浅 灰色 圆圈 表示 需求 量 ， 星 形 符号 表示 潜在 


设施 。 


图 1 案例 区 域 ZY( 左 ) 和 GY( 右 ) 
Figure 1 The study areas ZY (lest) and GY (right) 


首先 ， 基 于 两 个 地 理 区 域 数 据 构 造 两 个 测试 案例 。 每 个 区 域 中 ， 将 空间 单元 中 心 点 
位 置 作为 需求 点 位 置 和 候选 设施 位 置 ， 假 设 设施 i 服务 客户 j 的 费用 ci; = 1.0 * dij * dj» 
其 中 qj;j 为 设施 i 和 客户 j 间 的 欧 氏 距离 ， 单 位 为 公里 ，qj 为 客户 j 的 需求 量 。 设 施 建 设 
成 本 fi = (Ut 8i)si，5i 为 设施 服务 容量 ，s; 为 一 定 范围 内 的 随机 数 。 对 于 区 域 ZY， 取 
H=0.8、gi 为 区 间 [-0.1, 0.1] 内 的 一 个 随机 数 ， 对 于 区 域 ZY， 取 /=1.8、si 为 区 间 [-0.2, 0.2] 
内 的 一 个 随机 数 。 参 数 w 的 取 值 经 过 多 次 尝试 ， 使 SSCFLP 总 成 本 中 设施 成 本 比例 不 能 
太 高 ， 也 不 能 太 低 ， 处 于 30%~70% 之 间 。 

其 次 ， 基 于 ZY 和 GY 案例 数据 构造 测试 案例 集 。 针 对 每 个 案例 ， 扩 大 设施 服务 能 
力 为 1.2 倍 和 1.4 倍 ， 增 加 设施 成 本 为 1.1、1.2、1.3 和 1.4 倍 。 设 施 容 量 与 设施 成 本 组 
合 ， 每 个 区 域 上 共有 15 个 案例 。 案 例 基本 信息 见 表 1。 表 中 ,供需 比 等 于 所 有 候选 设施 
的 总 容量 与 总 需求 量 之 比 。 已 有 研究 表明 : 供需 比 是 案例 求解 难度 的 影响 因素 之 一 。 表 
中 ， 设 施 成 本 容量 比 等 于 候选 设施 的 总 成 本 与 总 容量 之 比 。 案 例 数 据 可 从 hittps://github. 
com/yfkong/unified 下 载 。 


表 1 测试 案例 集 


Table 1 The benchmark instances 


案例 I J| ”供需 设施 成 本 案例 I J 供需 ”设施 成 本 


比 容量 比 比 容量 比 
zyal 36 326 2.37 0.80 gyal 33 1276 1.61 1.83 
zya2 36 326 2.37 0.88 gya2 33 1276 1.61 2.02 
zya3 36 326 2.37 0.96 gya3 33 1276 1.61 2.20 
zya4 36 326 2.37 1.04 gya4 33 1276 1.61 2.38 
zya5 36 326 2.37 1.12 gya5 33 1276 1.61 2.56 
zybl 36 326 2.84 0.67 gybl 33 1276 1.93 1.53 
zyb2 36 326 2.84 0.74 gyb2 33 1276 1.93 1.68 
zyb3 36 326 2.84 0.80 gyb3 33 1276 1.93 1.82 
zyb4 36 326 2.84 0.87 gyb4 33 1276 1.93 1.98 
zyb5 36 326 2.84 0.94 gyb5 33 1276 1.93 2.14 
zycl 36 326 3.32 0.57 gycl 33 1276 2.26 1.31 
Zyc2 36 326 3.32 0.63 gyc2 33 1276 2.26 1.44 
zyc3 36 326 3.32 0.69 gyc3 33 1276 2.26 1.57 
zyc4 36 326 3.32 0.75 gyc4 33 1276 2.26 1.70 
zyc5 36 326 3.32 0.80 gyc5 33 1276 2.26 1.83 


3.2 案例 测试 


本 文 数学 启发 算法 使 用 Python 程序 设计 语言 编程 实现 。 算 法 步 又 如 下 : 读数 据 ,， 使 
用 ILH 方 法 获得 初始 解 ， 和 迭代 执行 邻 域 选 择 、 子 问题 求解 和 更 新 当前 解 ， 直 到 满足 终止 
条 件 。 子 问题 求解 环节 ， 使 用 PuLP(https://github.com/coin-or/pulp) 线 性 规划 建 模 工具 完 
成 SSCFLP 或 SSCKFLP 子 问题 模型 建构 ， 再 调用 IBM CPLEX 12.6 优化 器 完成 模型 求 
解 。 为 提升 计算 速度 ， 算 法 在 PyPy7 (https:Wwww.pypy.org) 环 境 中 运行 。 实 验 计 算 环 境 
为 : HP 桌面 计算 机 ， 配 置 Intel Core I7-6700 CPU 3.40-GHz 和 8GB 内 存 ，Windows 10 
操作 系统 。 本 文 算 法 可 从 网 页 https://github.com/yfkong/unified 下 载 ， 按 要 求 安装 软件 即 
可 运行 。 


= 


使 用 案例 数据 完成 4 个 问题 的 测试 。 所 有 案例 进行 测试 SSCFLP. CFLSAP 问题 求 
解 ; 针对 zyal 和 gyal 案例 ,测试 SSCKFLP、CKFLSAP 问题 求解 ， 案 例 zyal 设置 设施 
数量 K=13-22， 案 例 gyal 设置 设施 数量 K=18-30. 

使 用 两 种 方法 进行 问题 求解 : 精确 方法 和 本 文 数学 启发 方法 。 前 者 根据 本 文 模型 构 
造 案例 模型 ， 使 用 CPLEX 12.6 优化 嚣 求解。 记录 CPLEX 计算 结果 中 的 案例 目标 值 下 
界 、 目 标 值 上 界 及 计算 时 间 。 数 学 启发 算法 设置 参数 mloops 为 100。 一 般 来 说 ， 该 参数 
数值 越 大 ， 所 需 计算 时 间 也 越 长 。 记 录 计 算 结果 的 目标 值 、 目 标 值 与 下 界 差距 和 计算 时 
间 。 


3.3 计算 结果 


计算 结果 分 别 见 表 2、 表 3、 表 4 和 表 5。 表 中 ，LB 为 案例 目标 值 下 界 ，UB 为 案例 
目标 值 上 界 ，Gap 为 LB 和 UB 的 相对 差异 。 若 两 者 相等 表明 (Gap=0)， 表 明 CPLEX 找 
到 了 案例 最 优 解 。 表 中 最 右 三 栏 分 别 为 数学 启发 算法 获得 的 案例 目标 值 、 目 标 值 与 下 界 
相对 差异 Gap， 以 及 计算 时 间 。 


表 2 SSCFLP 计算 结果 
Table 2 Detailed SSCFLP solutions 

案例 CPLEX 数学 启发 算法 

LB UB Gap/% 时间 As 目标 值 Gap/% ”时 间 /s 
zyal 4428.23 4428.23 0.00 877.53 4428.23 0.00 62.23 
zya2 4722.23 4722.23 0.00 780.86 4722.23 0.00 72.15 
zya3 5025.23 5025.23 0.00 988.81 5025.57 0.01 58.55 
zya4 5327.60 5327.60 0.00 2042.05 5329.28 0.03 54.16 
zya5 5626.94 5626.94 0.00 1458.30 5628.48 0.03 88.73 
zybl 3999.99 3999.99 0.00 6024.45 4006.75 0.17 130.50 
zyb2 4239.99 4239.99 0.00 3368.16 4239.99 0.00 162.80 
zyb3 4488.99 4488.99 0.00 7106.34 4488.99 0.00 136.51 
zyb4 4736.99 4736.99 0.00 2676.81 4736.99 0.00 86.87 
zyb5 4982.99 4982.99 0.00 2662.05 4982.99 0.00 78.30 
zycl 3732.85 3732.85 0.00 227.04 3732.85 0.00 105.13 
zyc2 3940.85 3940.85 0.00 237.91 3940.85 0.00 89.73 
zyc3 4158.85 4158.85 0.00 457.83 4158.85 0.00 94.61 
zyc4 4374.85 4374.85 0.00 293.48 4374.85 0.00 96.94 
zyc5 4589.85 4589.85 0.00 758.59 4590.17 0.01 118.67 
gyal 3498793.22 3499012.55 0.01 7200.00 3499012.55 0.01 183.54 
gya2 3657078.11 3657323.63 0.01 7200.00 3657364.63 0.01 270.68 
gya3 3812527.81 3812605.72 0.00 7200.00 3812613.41 0.00 207.99 
gya4 3966893.30 3966921.73 0.00 7200.00 3966941.43 0.00 221.26 
gya5 4119130.39 4119543.12 0.01 7200.00 4119523.42 0.01 204.48 
gybl 3433138.33 3433223.99 0.00 7200.00 3433223.99 0.00 122.99 
gyb2 3573531.57 3573616.42 0.00 7200.00 3573616.42 0.00 151.29 
gyb3 3712631.49 3712713.37 0.00 7200.00 3717341.49 0.13 184.42 
gyb4 3847355.36 3847386.80 0.00 7200.00 3847386.80 0.00 194.18 
gyb5 3977990.12 3977990.12 0.00 2681.64 3978214.70 0.01 190.17 
gycl 3351292.23 3351292.23 0.00 587.87 3351292.23 0.00 203.33 
gyc2 3475209.95 3475209.95 0.00 863.11 3475209.95 0.00 144.90 
gyc3 3599127.67  3599127.67 0.00 1157.06  3599127.67 0.00 157.37 
gyc4 3723045.39 3723045.39 0.00 137.38  3724065.87 0.03 281.80 
gyc5 3840616.19 3840616.19 0.00 5198.12  3840702.03 0.00 313.58 


K 2 中 SSCFLP 计算 结果 表明 : CPLEX 能 够 求解 本 文案 例 ， 获 得 最 优 解 或 高 质量 
近似 解 (CPLEX MIPGap<0.01%)， 但 计算 时 间 差 异 较 大 ， a 2 a 最 慢 2 小 时 未 发 


现 最 优 解 。 数 学 局 发 算法 计算 速度 很 高 ， 各 案例 需要 1~6 分 钟 ， 计 算 结果 质量 很 高 ，12 
案例 获得 最 优 解 ， 其 余 案例 除 zybl 和 gyb3 外 ， 目 标 值 Gap 小 于 0.01%. 

K 3 CFLSAP 计算 结果 表明 : CPLEX 求解 本 文案 例 ， 均 获得 满意 的 结果 ， 平 均 Gap 
{EN 0.20%, 5 个 案例 在 2 小 时 时 间 内 获得 最 优 解 ; 与 表 2 相 比 , 增加 空间 连续 约束 后 ， 
CPLEX 计算 时 间 显 著 增加 。 数 学 启发 算法 需要 1-5 分 钟 时 间 求 解 各 案例 ， 平 均 Gap 值 
为 0.22%， 发 现 2 个 案例 最 优 解 ，12 个 案例 求解 质量 高 于 CPLEX 结果 。 

K 3 CFLSAP 计算 结果 
Table 3 Detailed CFLSAP solutions 


CPLEX MH 

案例 LB UB Gap/% ”时 间 /s 目标 值 ”Gap/% ”时 间 /s 

zyal 4435.21 4446.32 0.25 7200.00 4443.21 0.18 102.00 
zya2 4727.23 4740.69 0.28 7200.00 4740.18 0.27 211.85 
zya3 5031.18 5045.69 0.29 7200.00 5052.03 0.41 149.08 
zya4 5335.83 5344.66 0.17 7200.00 5365.55 0.56 62.45 
zya5 5634.87 5652.69 0.32 7200.00 5646.41 0.20 79.93 
zybl 4014.58 4014.58 0.00 3034.54 4015.33 0.02 92.11 
zyb2 4258.35 4263.26 0.12 7200.00 4263.08 0.11 233.67 
zyb3 4502.29 4519.58 0.38 7200.00 4528.97 0.59 111.21 
zyb4 4751.57 4778.24 0.56 7200.00 4779.08 0.58 139.89 
zyb5 4994.83 5032.58 0.75 7200.00 5030.37 0.71 173.50 
zycl 3736.16 3736.16 0.00 5988.96 3736.44 0.01 76.10 
zyc2 3941.78 3944.28 0.06 7200.00 3945.2 0.09 74.92 
zyc3 4161.52 4162.16 0.02 7200.00 4162.44 0.02 63.36 
zyc4 4372.17 4378.16 0.14 7200.00 4378.44 0.14 98.72 
zyc5 4593.16 4593.16 0.00 6698.44 4594.31 0.02 65.02 
gyal 3494350 3499875 0.16 7200.00 3499833 0.16 137.27 
gya2 3656614 3658269 0.05 7200.00 3658247 0.04 286.00 
gya3 3811904 3814304 0.06 7200.00 3815345 0.09 102.43 
gya4 3950935 3968368 0.44 7200.00 3968574 0.45 141.11 
gya5 4111963 4121919 0.24 7200.00 4120543 0.21 104.26 
gybl 3433408 3433408 0.00 5147.80 3433408 0.00 113.19 
gyb2 3573800 3573800 0.00 6247.53 3573800 0.00 135.53 
gyb3 3712195 3713060 0.02 7200.00 3712897 0.02 156.98 
gyb4 3844879 3852609 0.21 7200.00 3865510 0.54 130.93 
gyb5 3975431 3991801 0.41 7200.00 3978666 0.08 223.14 
gycl 3353103 3355378 0.07 7200.00 3357712 0.14 110.18 
gyc2 3475415 3482023 0.19 7200.00 3482121 0.19 117.73 
gyc3 3599353 3612109 0.35 7200.00 3605828 0.18 129.13 
gyc4 3723704 3730328 0.17 7200.00 3731717 0.22 239.15 
gyc5 3839448 3852543 0.34 7200.00 3850299 0.28 253.08 


K 4 中 SSCKFLP 计算 结果 表明 : 增加 设施 数量 约束 后 ，CPLEX 计算 时 间 大 幅 降 
低 ，20 个 案例 能 够 获得 最 优 解 ， 其 余 3 个 案例 获得 高 质量 解 (Gap<0.02%)，18 个 案例 仅 
需要 2 分 钟 计算 时 间 。 数 学 启发 算法 求解 质量 高 ， 发 现 11 个 案例 的 最 优 解 ，1 个 案例 目 
标 值 优 于 CPLEX 结果 ，Gap 均值 为 0.00%， 最 大 0.04%。 数 学 启发 算法 计算 时 间 处 于 
1~6 分 钟 。 


K 4 SSCKFLP 计算 结果 
Table 4 Detailed SSCKFLP solutions 


CPLEX 数学 启发 算法 
案例 K LB UB Gap/% ”时 间 /s 目标 值  Gap/% 时间 /s 
zyal 13 4846.88 4846.88 0.00 140.14 4848.83 0.04 101.15 
zyal 14 4740.14 4740.14 0.00 84.52 4740.42 0.01 107.92 
zyal 15 4644.42 4644.42 0.00 26.14 4644.47 0.00 88.33 
zyal 16 4549.83 4549.83 0.00 613.82 4549.83 0.00 193.39 
zyal 17 4479.05 4479.05 0.00 7.91 4479.05 0.00 75.13 
zyal 18 4459.38 4459.38 0.00 56.88 4459.38 0.00 56.50 
zyal 19 4428.23 4428.23 0.00 585.82 4428.57 0.01 62.06 
zyal 20 4441.19 4441.19 0.00 19.77 4441.19 0.00 74.95 
zyal 21 4490.37 4490.37 0.00 4.83 4490.37 0.00 34.00 
zyal 22 4617.61 4617.61 0.00 73.34 4617.88 0.01 44.80 
gyal 18 3516721.82 3517274.15 0.02 7200.00 3517195.85 0.01 329.68 
gyal 19 3501647.12 3501975.32 0.01 7200.00 3501975.32 0.01 203.00 
gyal 20 3498731.49 3499012.55 0.01 7200.00 3499012.55 0.01 140.92 
gyal 21 3499745.27 3499745.27 0.00 122.78 3499745.27 0.00 126.22 
gyal 22 3509618.91 3509618.91 0.00 63.27 3509618.91 0.00 87.52 
gyal 23 3531467.24 3531467.24 0.00 50.42 3531467.24 0.00 76.24 
gyal 24 355413439 355413439 0.00 40.91 355413439 0.00 88.87 
gyal 25 3584373.45 3584373.45 0.00 15.19 3584373.45 0.00 56.50 
gyal 26 3621833.78 3621833.78 0.00 14.58 3621840.41 0.00 44.42 
gyal 27 3663089.92 3663089.92 0.00 22.06 3663089.92 0.00 46.83 
gyal 28 3715849.01 3715849.01 0.00 18.59 3715872.84 0.00 47.81 
gyal 29 3772555.65 3772555.65 0.00 19.95 377257948 0.00 51.00 
gyal 30 3830365.47 3830365.47 0.00 14.66 3830366.59 0.00 51.07 


o 案例 计算 结果 表明 : 增加 空间 连续 约束 后 ，CPLEX 仍 能 够 高 质量 求 


解 本 文案 例 ， 但 计算 时 间 大 幅 增 加 ，16 个 案例 获得 最 优 解 ， 仅 2 个 案例 的 Gap 值 大 于 
0.10%. ay 启发 算法 能 够 快速 获得 高 质量 解 ， 共 发 现 6 个 案例 的 最 优 解 ，4 个 案例 优 于 
CPLEX 计算 结果 
K 5 CKFLSAP 计算 结果 
Table 5 Detailed CKFLSAP solutions 
案 = CPLEX 数学 启发 算法 
例 LB UB Gap/% ”时间 /s 目标 值 ” Gap/% 时间/s 
zyal 13 4855.13 4855.13 0.00 4816.65 4870.93 0.33 67.62 
zyal 14 4745.68 4745.68 0.00 3474.80 4760.25 0.31 138.40 
zyal 15 4650.36 4650.36 0.00 2200.54 4652.13 0.04 84.74 
zyal 16 4562.21 4562.59 0.01 7200.00 4570.24 0.18 80.41 
zyal 17 4486.15 4486.15 0.00 2200.07 4486.15 0.00 65.02 
zyal 18 4469.29 4469.29 0.00 5032.55 4475.10 0.13 67.40 
zyal 19 4438.07 4441.13 0.07 7200.00 4449.64 0.26 51.51 
zyal 20 4446.32 4446.32 0.00 2078.38 4447.13 0.02 59.24 
zyal 21 4500.53 4500.53 0.00 394.76 4506.77 0.14 24.85 
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zyal 22 4619.83 4620.24 0.01 7200.00 4632.94 0.28 28.21 
gyal 18 3513204.32 3522650.84 0.27 7200.00 3518302.79 0.15 53.45 
gyal 19 3501266.20 3509894.00 0.25 7200.00 3502926.68 0.05 115.24 
gyal 20 3498654.87 3500110.57 0.04 7200.00 3499826.97 0.03 131.71 
gyal 21 3500566.31 3500566.31 0.00 4194.36 3500566.31 0.00 95.69 
gyal 22 3510439.96 3510439.96 0.00 4944.07 3510439.96 0.00 67.28 
gyal 23 3532049.19 3532374.32 0.01 7200.00 3532288.28 0.01 82.74 
gyal 24 3554955.43 3554955.43 0.00 1935.15 3554955.43 0.00 64.27 
gyal 25 358519449 3585194.49 0.00 803.02 3585251.36 0.00 55.72 
gyal 26 3622654.82 3622654.82 0.00 1113.47 3622654.82 0.00 54.34 
gyal 27 3663910.96 3663910.96 0.00 343.29 3663910.96 0.00 54.92 
gyal 28 3716663.42 3716663.42 0.00 591.89 3716714.23 0.00 43.74 
gyal 29 3773370.06 3773370.06 0.00 405.58 3773371.19 0.00 50.81 
gyal 30 3831365.77 3831365.77 0.00 120.78 3831393.55 0.00 45.43 
3.4 比较 分 析 


比较 SSCFLP 和 CFLSAP 计算 结果 ， 有 具有 以 下 特征 : (1)SSCFLP 增加 设施 服务 区 空 
间 连 续 约束 ， 目 标 值 略 有 增加 ，ZY 案例 增加 0.09%~1.00%， 平 均 0.39%; GY 案例 增加 
0.01%~ 0.36%， 平均 0.13%。(2)SSCFLP 增加 空间 连续 约束 后 ，CPLEX 求解 时 间 显 著 增 
加 , 求解 质量 也 有 下 降 , 目标 值 与 下 界 差异 MIPGap 从 0.00% 增 加 到 0.20%. (3) SSCFLP 
增加 空间 连续 约束 ， 对 数学 启发 算法 计算 时 间 无 明显 影响 。(4) SSCFLP 增加 空间 连续 约 
束 对 设施 区 位 和 服务 区 有 显著 影响 。 本 文 测试 案例 中 ， 部 分 案例 无 论 是 否 考 虑 服务 区 空 
间 连 续 ， 选 择 设 施 区 位 完全 相同 ， 但 另 一 部 分 案例 ，SSCFLP 和 CFLSAP 区 位 选择 有 显 
著 差 异 。 例如， 案例 zya4， 图 2 为 SSCFLP 和 CFLSAP 计算 结果 示意 图 。 两 者 目标 值 差 
异 0.39%， 设 施 数量 均 为 19， 但 3 个 设施 区 位 有 差异 ; 前 者 4 个 设施 的 服务 区 不 连续 ， 
后 者 保证 所 有 服务 区 空间 连续 。 


图 2 案例 zya4 计算 结果 ( 左 为 SSCFLP， 右 为 CFLSAP) 
Figure 2 SSCFLP (left) and CFLSAP (right) solutions on the instance zya4 
设施 数量 约束 对 于 案例 目标 值 影响 显著 。 首 先 ， 无 论 SSCFLP 或 CFLSAP 问题 ， 增 
加 设施 数量 约束 后 ， 精 确 算法 的 求解 时 间 显著 下 降 。 可 能 的 原因 是 : 增加 设施 数量 约束 
后 ， 可 行 的 候选 设施 组 合 数量 大 幅 降 低 ， 大 幅 缩 小 了 可 行 解 空间 。CPLEX 计算 结果 表 
明 ， 增 加 设施 数量 约束 后 ，SSCKFLP HE SSCFLP 更 容易 获得 最 优 解 ，CKFLSAP 比 
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CFLSAP 容易 获得 最 优 解 。 其 次 ， SSCFLP 或 CFLSAP 存在 最 优 设施 数量 , 增加 设施 数 
量 约束 ， 设 施 数量 偏离 原 问题 最 优 设施 数量 越 多 ， 目 标 增加 值 越 多 。 图 3 为 SSCKFLP 
目标 值 随 设施 数量 变化 曲线 ， 上 图 为 案例 zyal， 下 图 为 案例 gyal， 两 个 案例 SSCFLP 最 
好 解 中 ， 设 施 数量 分 别 是 19 和 20。 
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图 3 SSCKFLP 目标 值 ( 纵 轴 ) 随 设施 数量 ( 横 轴 ) 变 化 曲线 
Figure 3 SSCKFLP objectives on instances zyal (upper) and gyal (lower) versus the number 
of facilities 

在 需求 分 布 确定 、 候 选 设施 位 置 确定 的 情况 下 ， 设 施 容 量 和 固定 成 本 对 区 位 选择 有 
一 定 的 影响 。 本 文案 例 中 ， 按 同一 比例 调整 设施 容量 ， 设 施 区 位 选择 通常 会 发 生变 化 。 
而 同一 比例 调整 设施 成 本 ， 部 分 案例 的 区 位 选择 仍 保持 不 变 。 比 较 发 现 ， 以 下 各 组 案例 
的 设施 区 位 选择 完全 相同 : zyal、zya2 和 zya3, zybl, zyb2, zyb 3、zyb 4 和 zyb 5, 
Zzycl、zyc2、zyc3、zyc4 和 zyc5, gya2 和 gya3, gya4 和 gyaS, gyb2 和 gyb3, gycl. gyc2 
和 gyc3。 


4 结论 

本 文 对 经 典 的 SSCFLP 进行 了 扩充 ， 即 增加 设施 服务 区 空间 连续 约束 和 设施 数量 约 
R, 衍生 出 SSCKFLP、CFLSAP 和 CKFLSAP 这 3 个 变种 问题 。 利 用 设施 区 位 ， 借 鉴 网 
络 流 模 型 (Duque et al. 2011)， 作 者 构建 了 顾及 空间 连续 的 CFLSAP 和 CKFLSAP 模型 。 


使 用 流 模型 构建 p-regions 区 划 问 题 模型 , 模型 精确 求解 仅 适用 于 规模 很 小 的 案例 ; 而 本 
X CFLSAP 和 CKFLSAP 模型 相对 容易 求解 , 原因 在 于 SSCFLP 选择 的 设施 区 位 可 以 直 
接 用 于 建立 流 模型 。 实 验 表 明 : CPLEX 优化 器 能 够 求解 1276 个 单元 和 33 个 设施 的 案 
例 ， 获 得 最 优 解 或 高 质量 解 (MIPGap< 0.44%). SSCFLP 应 用 实践 中 ， 可 能 要 求 设 施 服务 
区 连续 ， 也 可 能 需要 指定 设施 数量 ， 本 文 SSCFLP 扩展 模型 符合 这 些 应 用 需求 。 

本 文 设计 了 数学 启发 算法 求解 SSCFLP、SSCKFLP、CFLSAP 和 CKFLSAP。 数 学 
启发 算法 本 质 上 是 超大 规模 邻 域 搜索 ， 即 从 初始 解 作 开始 ， 使 用 数学 模型 迭代 地 改进 当 
前 解 中 的 超大 国 模 邻 域 ， 直 到 当前 解难 以 更 新 。 案 例 测试 表明 ， 数 学 启发 算法 求解 
SSCFLP、SSCKFLP、CFLSAP 和 CKFLSAP， 目 标 值 与 其 下 界 的 差异 分 别 为 0.01%、 
0.00%、0.22% 和 0.08%， 求 解 质量 与 CPLEX 计算 结果 无 明显 差异 ， 但 计算 效率 大 幅 提 
升 。 案 例 测试 验证 了 本 文 数学 启发 算法 的 有 效 性 。 

本 文案 例 研究 也 发 现 : SSCFLP 或 SSCKFLP 增加 设施 服务 区 空间 连续 约束 后 ， 目 
标 值 增加 幅度 不 大 ， 但 最 优 设施 区 位 可 能 发 生变 化 。SSCFLP 或 CFLSAP 增加 设施 数量 
约束 对 于 目标 值 影响 显著 ， 设 施 数量 偏离 原 问 题 最 优 设施 数量 越 多 ， 目 标 增加 值 越 多 。 
在 需求 分 布 确定 、 候 选 设施 位 置 确定 的 情况 下 ， 设 施 容 量 和 固定 成 本 对 区 位 选择 有 一 定 
的 影响 。 
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A matheuristic algorithm for the single-source capacitated facility 
location problem and its variants 

Abstract: This article proposes a matheuristic algorithm for the single-source facility location 
problem (SSCFLP) and its variants: SSCFLP with K facilities (SSCKFLP), SSCFLP with 
connective service areas (CFLSAP), and SSCFLP with K facilities and connective service areas 
(CKFLSAP). The algorithm starts from an initial solution, and then iteratively improves the 
solution by searching large-scale neighborhood of current solution. The neighborhood is defined 
by determining a subset of candidate facilities and a subset of customers: (1) randomly select a 
customer; (2) select Q nearest facilities and their customers from the current solution; (3) select 
nearest candidate facilities of all the customers; and (4) randomly drop some candidate facilities 
if too many facilities are selected in previous step. The size of neighborhood is critical to the 
performance of the algorithm: it is hard to solve an extra-large neighborhood and it is difficult 
to find a better solution in a small neighborhood. The value of Q is suggested according to the 
number of facilities in current solution. The current solution might be improved by finishing the 
following steps: (1) formulate a sub-problem model using the selecting facilities and customers; 
(2) solve the model and update the current solution using sub-problem solution; and (3) for 
CFLSAP or CKFLSAP, repair the non-connective service areas, and improve solution with local 
search operators. Two set of instances were generated to test the algorithm. Experimentation 
shows that the instances of SSCFLP and its variant problems can be solved by the proposed 
matheuristic algorithm effectively and efficiently. The solutions found by the proposed 
algorithm approximate optimal solutions or the lower bounds with average gaps of 0.01% for 
SSCELP, 0.22% for CFLSAP, 0.00% for SSCKFLP, and 0.08% for CKFLSAP. Solution results 
show that the solution objective would be slightly increased by adding the contiguity constraints 
on SSCFLP or SSCKFLP. The optimal facility locations of SSCFLP/SSCKFLP might be 
different from those of CFLSAP/CKFLSAP. 

Keywords: single-source capacitated facility location problem; connective service area; 
matheuristic algorithm; Large-scale neighborhood. 
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